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ABSTRACT 

We begin with a review of the antiferromagnetic spin 1/2 Heisen- 
berg chain. In particular, we show that the model has particle-like 
excitations with spin 1/2, and we compute the exact bulk S matrix. 
We then review our recent work which generalizes these results. We 
first consider an integrable alternating spin 1/2 - spin 1 chain. In 
addition to having excitations with spin 1/2, this model also has exci- 
tations with spin 0. We compute the bulk S matrix, which has some 
unusual features. We then consider the open antiferromagnetic spin 
1/2 Heisenberg chain with boundary magnetic fields. We give a direct 
calculation of the boundary S matrix. (Talk presented at the confer- 
ence on Statistical Mechanics and Quantum Field Theory at USC, 16 
- 21 May 1994) 



1. Introduction 

The investigation of integrable quantum spin chains was initiated by Bethe 1 
with the classic paper on the closed * spin 1/2 Heisenberg chain. Other examples 
of integrable quantum spin chains include the open spin 1/2 chain 2 " 4 , the spin 1 
chain 5 ' 6 , the spin 1/2 chain with a spin 1 impurity 7 , and the alternating spin 1/2 - 
spin 1 chain 8 . 

There are several motivations for studying integrable quantum spin chains. 
First, these are many-body quantum mechanical models for which exact results 
can be computed. Also, these models typically have 9 ' 10 a regime with a nontrivial 

* One-dimensional quantum spin chains, like strings, come in two topologies: 
closed (periodic boundary conditions) and open. 



1 



antiferromagnetic vacuum and novel excitations ("spinons"). In the continuum 
limit, these excitations are described by 1 + 1-dimensional integrable relativistic 
quantum field theory. (See, e.g., Refs. 11-14.) Last but not least, such models have 
applications in statistical mechanics and condensed matter physics 15 and perhaps 
also in string theory 16 . 

In this talk we review our recent work on both bulk and boundary S matrices 
for the excitations of integrable quantum spin chains. Such S matrices provide 
valuable information about long-distance physics and boundary phenomena of the 
models. (See, e.g., Refs. 17, 18 and references therein.) 

The outline of this talk is as follows. We begin with a brief review of the closed 
spin 1/2 Heisenberg chain in the antiferromagnetic regime, with emphasis on the 
physical properties which emerge from the Bethe Ansatz solution. In particular, 
following Faddeev and Takhtajan 9 , we outline the argument that the model has 
particle-like excitations with spin 1/2. These excitations interact, and we explain 
how the exact S matrix can be computed. 

In the remainder of the talk, we generalize these results in two different di- 
rections. In Section 3 we consider the alternating spin 1/2 - spin 1 chain in the 
antiferromagnetic regime. Following Ref. 19, we show that in addition to having 
excitations with spin 1/2 (as in the Heisenberg chain), this model also has excita- 
tions with spin 0. We compute the S matrix, which has some unusual features. In 
Section 4, we consider the open antiferromagnetic spin 1/2 Heisenberg chain with 
boundary magnetic fields. Following Ref. 20, we give a direct calculation of the 
boundary S matrix. This is the first first-principles calculation of a boundary S 
matrix corresponding to an interacting relativistic field theory. Our result agrees 
with the boundary S matrix for the boundary sine-Gordon model with (3 2 — > 8% 
and with "fixed" boundary conditions 18 ' 21 . 



2. Closed Spin 1/2 Chain 

The Hamiltonian of the closed antiferromagnetic isotropic spin 1/2 Heisenberg 
chain is given by 

1 N 

Ti = -^2(a n ■ a n+ i - 1) , a N+1 =a 1 , (2.1) 

n=l 

where a are the usual Pauli spin matrices. We assume that the number of spins, 
N, is even. The Hamiltonian commutes with the "momentum" operator P (defined 
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such that e lP is the one-site shift operator), as well as with the sit (2) generators 
S = \ J2n=i ^n- The so-called Bethe Ansatz states are the simultaneous eigenstates 
of Ti, P, S 2 and S z which are highest weights of su(2) (i.e., with corresponding eigen- 
values S = S z > 0). These states have been determined by both the coordinate 1 
and algebraic 22 Bethe Ansatz methods.* In the latter approach, one constructs 
certain creation and destruction operators, B(X) and C(A), respectively; and the 
Bethe Ansatz states are given by 

B(X 1 ) B(X 2 )---B(X M ) u+, (2.2) 

where u + is the ferromagnetic vacuum state with all spins up, 

C(A)w+=0, (2.3) 

and {A Q } satisfy the Bethe Ansatz (BA) equations 

The corresponding eigenvalues are given by 

M 



1 1 
E -~o l^YTT 1 



M 




(2.5) 



For the ferromagnetic spin chain with Hamiltonian — Ti, the ground state has 
all spins aligned, and evidently corresponds to M = 0. For the antiferromagnetic 
spin chain with Hamiltonian +Ti (2.1), the identification of the ground state and 
the lowest-lying excited states is as follows: 



* The remaining states are obtained by acting on the Bethe Ansatz states with 
the spin lowering operator S~ . 
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2.1 Ground state 

For the ground state, one can argue that M = ^ and that the roots {Ai , • • • , Am} 
are all distinct and real. This solution of the BA equations corresponds to a filled 
Fermi sea (i.e., no holes). See Fig. 1. 

For N — > oo, the set of A's becomes dense on the real line, and is described by 
the density a vac (\) which is given by* 

(This result is obtained by solving the linear integral equation for the root density 
which follows from the BA equations.) Making in Eq. (2.5) the following replace- 
ment of sums by integrals 

A? £()->/ (K-c(A)dA, (2-7) 

a=l J -°° 

one concludes that the ground state has the following quantum numbers: 

E = E = —N log 2 , P = P Q = ATtt/2, S = 0. (2.8) 

In particular, the ground state is a spin singlet, as one would expect for an antifer- 
romagnet. 

2.2 Excitations 

The excited states above the ground state consist of an even number of particle- 
like excitations, which are now known as "spinons". (Faddeev-Takhtajan called 
them "kinks".) Therefore, the lowest-lying excited states have two spinons. One 
can argue that there are four such states: the triplet (S = 1) states, and the singlet 
(S = 0) state. (The total number of states with v spinons is equal to 2 V .) The 
fact that the excited states with two spinons have S = 1 and S = implies the 
important result that a spinon has spin 1/2. 

The triplet state with S = S z = 1 is described by only real roots {Ai , • • • Am} 
(as is the ground state), but with M = -y — 1. This solution of the BA equations 
corresponds to a Fermi sea with two holes. The hole rapidities are labeled Ai, A 2 . 
See Fig. 2. 

* The Fermi points are at ±oo. Had we introduced a bulk magnetic field, the 
Fermi points would be at A = ±A, with A finite. 
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The singlet state corresponds to a Fermi sea with two holes (with rapidities A\ 
and A 2 ) as well as a "2 - string", which is a set of conjugate roots A ± §, with Ao 
(the "center" of the 2 - string) real. For the singlet state, the BA equations further 
constrain the center to be given by 

A = ^(Ai + A 2 ). (2.9) 

All of these features of the singlet-state solution can be seen in Fig. 3. 

For both the triplet and singlet states, one can show that the density cr(A) of 
real roots and holes is given by an expression of the form 

where r(A) is a correction of order 1 to the ground state density (2.6). Heuristically, 
this correction corresponds to a "polarization" of the Fermi sea due to the holes 
and the 2-string. Explicitly, 



^triplet(^) — J(X — A Q ) , 

a=l 

TsingletW = TtripletW ~ «l(A - A ) , ( 2 -ll) 

where 

i r°° p~ M i rl ( \ -\- i- \ 

J(x) = rJ_^^ x ih^' -w = dH^tj- (2 ' 12) 

It follows from Eq. (2.5) that for both the triplet and singlet states, the energy and 
momentum are given by 

£ = £ +£(Ai)+£(A 2 ), 

P = P +p(Ai)+p(A 2 ), (2.13) 

where Eq and Pq are the energy and momentum of the ground state, and 

e(A) = - - , (2.14) 

v ' 2cosh7rA v ' 

p(A) = tan _1 sinh7rA - - . (2.15) 



5 



From the additivity property displayed by Eq. (2.13), we see that the spinons indeed 
are particle-like excitations, with energy e(X) and momentum p(X). The energy- 
momentum dispersion relation is 

71" 

e = — — sinp . (2-16) 

Note that the spinons are gapless (e(X) — > for A — > ±oo). 
2.3 S matrix 

The S matrix for the scattering of spinons can be calculated exactly. Here we 
follow the Korepin-Andrei-Destri 23 ' 24 method. An important observation is that 
for a state of two spinons with rapidities Ai and A 2 , the momentum p(Ai) satisfies 
the quantization condition 

e *p(>i)JV ^(A x - A 2 ) = 1 , (2.17) 

where R is the 2-particle S matrix (acting in the tensor product space C 2 <E> C 2 ), 
and N is the number of spins in the chain. Let e 1 ^ be an eigenvalue of R. Then 
p(Xi) is related to the phase shift 4> by 

p(A 1 ) + l0=^m, (2.18) 



where m is an integer. 

On the other hand, one can show that 

2tt [~ Xl 2tt ~ 

p(Ai) + — J r(X) dX + const = —J lt (2.19) 

where r(A) is the function appearing in Eq. (2.10), and J\ is an integer or half-odd 
integer. Comparing Eqs. (2.18) and (2.19), we conclude that the phase shift <j> is 
given by 

/■Ai 

<p = 2n r(A) dX + const . (2.20) 



Using the explicit expressions for r(A) for the triplet and singlet states, we obtain 
(up to a rapidity-independent phase factor) 

triplet (A) = e^ut = iU+ 2 / 12 2^ , 

p n ' r(i-f )r(| + f ) 

Ssinglet{X) = e*'™'* = -^S triplet (X) , (2.21) 

where A = A\ — A 2 . 
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2.4 Further remarks 

It is useful to formulate the above result as a 4 x 4 matrix. Since R commutes 
with su(2), it is a linear combination of the identity matrix 1 and the permutation 
matrix V. Moreover, R has the eigenvalues (2.21). It follows that 



The matrix R satisfies unitarity and crossing, as well as the Yang-Baxter equation 



where _Ri 2 , -R13, and R23 are matrices acting in the tensor product space C 2 ®C 2 <g> 
C 2 , with R12 = R <g) 1, -R23 = 1 <8> .R, etc. (See, e.g., Refs. 17, 25 and references 
therein.) 

We have described here only the calculation of the 2-particle S matrix. In prin- 
ciple one can compute in similar fashion the multiparticle S matrix, and verify that 
the multiparticle S matrix is factorizable into a product of 2-particle S matrices. 

Finally, we briefly discuss the continuum limit of this model. The continuum 
quantum field theory is 26 the su(2) WZW model 27 of level k = 1. This is an su(2)- 
invariant CFT 28 with central charge c = 1. Indeed, the massless S matrix 13 of the 
latter model coincides with Eq. (2.21). 

The S matrix (2.21) can also be obtained by starting with an anisotropic spin 
chain with anisotropy parameter 77 and lattice spacing a, and then taking the con- 
tinuum limit a — > and the isotropic limit 77 — > while keeping a mass parameter 
m 2 oc a -2 exp(— 7r 2 /?7) fixed. (See Ref. 11.) Thus, this S matrix also describes 
a massive sw(2)-invariant integrable quantum field theory (namely 17 , the su(2)- 
invariant Thirring model, or the sine-Gordon model in the limit (3 2 — > Sir), which 
in the ultraviolet limit m — > reduces to the WZW model. 

3. Closed Alternating Spin 1/2 - Spin 1 Chain 

We now consider a system with a strictly alternating arrangement of 2N spins, 
with spins 1/2 at even sites and spins 1 at odd sites. That is, there are iV spins 
\&2 7 1<?4 » • • • > \&2N °f s P m 1/2 and iV spins s 1 , s 3 , • • • , sW-i of spin 1. The 
sit (2) -invariant Hamiltonian H is given by 8 ' 29 



R(X) = VR(X) 



r(i + f)r(i-f)(Ai-zP) 
r(i-f)r(i + f) (A-i) 



(2.22) 



R 12 (X - X') R 13 (X) R 23 (X') = R 23 (X') R 13 (X) R 12 (X - X') , 



(2.23) 




(3.1) 



n=l 
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+ (2<7 2 n • S 2n -1 + 1) [(1 + S 2 n-1 ' S 2n +l) (2<72n ' «2n+l + 1) + 2] | . 

Note that the Hamiltonian contains both nearest and next-to-nearest neighbor in- 
teractions. We assume periodic boundary conditions {fii n = &2n+2N and S2n+i = 
S2u+i+2n) and that N is even. 

It is not difficult to understand the origin of this Hamiltonian. Using R ma- 
trices* fl(3>3), Rih 1 ), Ri 1 ^) and i?^ 1 ' 1 ) as vertex weights, one can construct an 
integrable two-dimensional classical statistical mechanical vertex model as shown 
in Fig. 4. Note that both rows and columns alternate between spin 1/2 and spin 1, 
and that the lattice is invariant under rotation by tv/2. The logarithmic derivative 
of the (two-row to two-row) transfer matrix gives the above Hamiltonian. 

The Bethe Ansatz states have been determined in Ref. 8. The corresponding 
energy, momentum, and spin eigenvalues are given by 

M (\ 1 1 \ 

E = -J2 [2 \2 + 1 + \2 + 1 J + ^dependent of {A a } , 

P=-yio g ( Xa + ^. X - + l ) , (3.2) 

S = S>= 3 4-M, 
2 

where the variables A a satisfy the BA equations 



JV 



M 



A + |A a+ » = X x e + i M<f.(3.3) 
A Q - f A Q - 1 J a± \ a _ \p _ % 2 



0^a 



The momentum operator is defined such that e l2P is the two-site shift operator, 
and hence the factor 1/2 in Eq. (3.2). 

The ground state corresponds to two filled Fermi seas: a sea of 1-strings (i.e., 
real roots of the BA equations, as in the ground state of the spin 1/2 chain) and a 
sea of 2-strings. See Fig. 5. 

Holes in the sea of 2-strings are excitations with spin 1/2, just like the ex- 
citations of the spin 1/2 chain. However, for the alternating spin chain, there is 
the additional possibility of having holes in the sea of 1-strings. As shown in Ref. 



* We denote by R( Sl ' S2 - ) (A) the sw(2)-invariant R matrix acting on the tensor 
product space C 2si+1 <8> C 2s2+1 corresponding to spins s\ and s 2 - See Refs. 5, 30. 
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19, holes in the sea of 1-strings are excitations with spin 0. To the best of our 
knowledge, this is the first example of a magnetic chain with spin excitations. 

For both the spin 1/2 and spin excitations, the energy e(X) is given by 
(2.14), and the momentum is given p(X)/2, where p(X) is given by (2.15). The total 
number of states with v excitations is X]m=o^ 2m ? which corresponds to having an 
even number of each type of excitation. 

The S matrix can be computed (up to rapidity-independent phase factors) as 
before. The triplet and singlet S matrix elements for the scattering of two spin 1/2 
excitations coincide with the expressions given in Eq. (2.21). There is no scattering 
between two spin excitations (the S matrix element is S(X) = 1) and the S matrix 
element for the scattering of a spin 1/2 excitation and a spin excitation is 

S(A)=tcoth!(A+l) , (3.4) 

where A is the difference of the corresponding hole rapidities. Remarkably, the 
scalar-spinor scattering is nontrivial, yet the spinor-spinor scattering is the same as 
for the Heisenberg chain. 

An interesting open problem is to determine the continuum limit of this model. 
We know that the continuum quantum field theory must be some sw(2)-invariant 
CFT with 29 ' 31 central charge c = 2. 

We remark that for both the spin 1/2 chain and the alternating spin 1/2 - 
spin 1 chain, the ratio Cr/T (the specific heat at constant field divided by the 
temperature) has the property 

lim lim = lim lim — ^- . (3.5) 

The LHS can be evaluated by the method of Filyov, et a/. 32 while the RHS can be 
evaluated by the method of Johnson and McCoy 33 . 

For integrable isotropic spin s chains with s > 1/2 30 > 34 > 35 5 the property (3.5) is 
not satisfied. Indeed, the LHS is proportional to c = 3s/ (s + 1) (see Ref. 26), while 
the RHS is proportional to c = 1 (see Ref. 36). Moreover, there is a discrepancy 
between the results of Takhtajan 34 (see also Ref. 24) and Reshetikhin 37 for the 
two-body S matrix: 

STakhtajan 7^ $ Reshetikhin ■ (3-6) 

These facts strongly suggest that there are (at least) two continuous field theories 
in the (T , H ) = (0,0) limit of the spin s isotropic chain. The limit T = , H = + 
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corresponds to a c = 1 theory with Takhtajan's S matrix; and the limit H = , T = 
+ corresponds to a c = 3s/ (s + 1) theory with Reshetikhin's S matrix. 

4. Open Spin 1/2 Chain with Boundary Magnetic Fields 

We now consider the open antiferromagnetic isotropic spin 1/2 Heisenberg 
chain with boundary magnetic fields. The Hamiltonian is given by 



where the (real) parameters £± correspond to boundary magnetic fields. We assume 
that £± > 1/2 and that N is even. Since the spin chain is open, the Hamiltonian 
does not commute with the shift operator. Moreover, the boundary magnetic fields 
break the su(2) symmetry, and so the Hamiltonian commutes only with S z . 

The simultaneous eigenstates of Ti and S z have been determined by both the 
coordinate 3 and algebraic 4 Bethe Ansatz. In the latter approach, one constructs 
(in analogy with the closed spin chain) certain creation and destruction operators, 
£>(A) and C(A), respectively; and the eigenstates are given by 




(4.1) 



B(Ai) B{M)---B{\m) c^ + , 



(4.2) 



where u + is the ferromagnetic vacuum state with all spins up, 



C(A) u + =0, 



(4.3) 



and {A Q } satisfy the Bethe Ansatz (BA) equations 






(4.4) 



The corresponding energy and spin eigenvalues are given by 



S z 




(4.5) 



(4.6) 
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We require that the BA solutions correspond to independent BA states, and there- 
fore, we make the restriction 

Re(A a )>0. (4.7) 

(See, e.g., Refs. 2, 3, 20, 21, 38, 39.) 

As for the closed spin 1/2 chain, the ground state corresponds to a real Fermi 
sea, and the excitations are spinons with S z = ±1/2 and energy e(X). We assume 
that the 2-particle S matrix is the same as for the closed spin chain. The problem 
is to compute the boundary S matrix, which describes the interaction of a spinon 
with the end of the spin chain. However, it is instructive to first consider a similar 
but more elementary problem. 

4-1 Boundary S matrix: free particle 

As a warm-up exercise, we first compute the boundary £ matrix for a free non- 
relativistic particle of mass m (with Hamiltonian H = p 2 /2m) which is constrained 
to be on the positive half-line x > 0. Usually one demands that the wavefunction 
ip(x) vanish at x = 0. This is a sufficient, but by no means necessary, condition 

for the probability current j(x) = itfj(x)* d x tfj(x) to vanish at x = 0. We consider 
instead the more general (mixed Dirichlet-Neumann) boundary condition 

d 

cibix) + -rip(x) = at x = , (4.8) 
dx 

where c is a real parameter with dimension 1 /length. This boundary condition also 
implies the vanishing of the probability current at x = 0, and is compatible with the 
self-adjointness of the Hamiltonian. (This boundary condition has been shown 40 to 
be compatible with the integrability of the nonlinear Schrodinger equation on the 
positive half-line.) Assuming energy eigenfunctions of the plane- wave form 

iP P (x) = Ae ipx + Be~ ipx (4.9) 

(we set h = 1), we can use the boundary condition (4.8) to eliminate A in terms of 
B; and we immediately obtain 



ip p (x) = B 



p + ic 



-ipx _|_ f > e ipx 

K p-ic 



(4.10) 



We conclude that the boundary £ matrix is given by 



K(p) = P±^. (4.11) 

p — ic 
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We see that the boundary can give rise to a nontrivial boundary S matrix. The 
pole at p = ic implies the existence (for c > ) of a boundary bound state with 
energy E = — c 2 /2m. 

4-2 Boundary S matrix: open spin 1/2 chain 

For the open spin 1/2 chain, the quantization condition (2.17) is replaced by 21 

e i2p(Ai)Ar Ri2{ ~ Xi _ ~ X2) Ki(Ai ,£_) R 21 (\i + A 2 ) K^Xi ,£+) = 1 • (4.12) 

Here p(A) is defined by (2.15) (i.e, the expression for the momentum of a particle 
with rapidity A for the corresponding system with periodic boundary conditions), 
and K(\,£) is the boundary S matrix (acting in the space C 2 ). We use the same 
notation employed in Eq. (2.23); moreover, 



i?2l(A) = Pi2 Ri2(X) V12, 



(4.13) 



where V is the permutation matrix; and Ki, K 2 denote matrices acting in the space 
C 2 <g> C 2 , with K 1 =K®1,K 2 = 1®K. 

The R matrix is given in Eq. (2.22). This matrix has the following form 



R(X) = 



a(X) 








° \ 





6(A) 


c(A) 








c(A) 


6(A) 














a(X)J 



(4.14) 



with 

6(A) 



A 



X-i 



to(A) , c(A) 



X-i 



a(A) , a(A) 



r(i + f)r(i-f) 



The U(l) symmetry of the Hamiltonian's boundary terms implies that the 
boundary S matrix is of the form 



K(X,0=[[ ai Y ) 





/5(A,0 



(4.15) 



Our task is to explicitly determine the matrix elements a(X , £) and (3(X , £), which 
are the boundary scattering amplitudes for excitations with S z = +1/2 and S z = 
—1/2, respectively. We proceed by examining the two-particle excited states, which 
we classify by their S z eigenvalue. As for the closed spin 1/2 chain, there are 
four such states (S z = 1, S z = — 1, and two states with S z = 0). Since we need 
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to determine only two matrix elements, the system of four equations provided by 
the quantization condition (4.12) is overdetermined. The structure (4.14) of the R 
matrix suggests that there will be two simple relations corresponding to the diagonal 
elements of the R matrix. These relations will enable us to determine the matrix 
elements a(A,£) and /?(A,£). The other two relations should lead to identities. 

S z = 1 state 

For the S z = 1 state, the quantization condition (4.12) implies 

2p(\ 1 ) + -*M = -m, (4.16) 



with 

e** (1) = a (Ai - A 2 ) a(X 1 , £_) a(X 1 + A 2 ) a(X 1 , £+) • (4.17) 



As for the closed spin chain, the S z = 1 state is the Bethe Ansatz state consist- 
ing of two holes in the (real) Fermi sea. Using the BA equations, we can compute 20 
the function r(A), which is the sum of 1/N contributions to the density er(A) for 
this state.* For the open spin chain, the identity (2.19) is replaced by 

2n f^ 1 2n ~ 

2p(X 1 ) + — J r(X) dX + const = — Ji . (4.18) 

It follows that 

= 2tt f 1 r(A) dX + const . (4.19) 
Jo 

Using the explicit expressions for r(A) and a(A), we obtain the following result for 
a(X,£) (up to a rapidity-independent phase factor): 

n _ r(^ + j) r(f + i) r(^ + i(2^-i)) r(f + + 
r (t + l) r (^ + !) r(f + 1(2^-1)) r(^A + i(2e + i))- 

(4.20) 

S 2 = -1 state 

To determine the remaining element /3(A,£) of the boundary S 1 matrix, we 
consider the S z = —1 state. The quantization condition (4.12) implies 

2p(A 1 ) + l$(- 1 ) = ^m, (4.21) 



* In contrast to the closed-chain result (2.6), the ground-state density cr„ ac (A) for 
the open spin chain has corrections of order 1/N, and these corrections contribute 
to r(A). 
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with 

e^ (_1) = a(X 1 - A 2 ) P(X 1 , £_) a(Ai + A 2 ) P(X 1 , £+) . (4.22) 



The S z = — 1 state is most easily described within the BA approach by changing 
the pseudovacuum. Hence, instead of working with the states (4.2), we work now 
with 

C(Ai) C(X 2 )---C(X M ) oo- , (4.23) 
where u~ is the ferromagnetic vacuum state with all spins down, 

B(A)uT = 0. (4.24) 

Sklyanin has shown 4 that {A a } in Eq. (4.23) satisfy the same BA equations (4.4) as 
before, except for the replacement of £± by — £±. The energy eigenvalues are given 
by the same expression (4.5), and the S z eigenvalues are now given by 

N 

S z = M - — . (4.25) 

The S z = — 1 state now corresponds to the Bethe Ansatz state consisting of 
two holes in the Fermi sea. The calculation of the function r(A) is exactly the same 
as for the S z = 1 state, except that we must track the change £± — > — £±. We find 
that j3{X , £) is given by 

/3(A,0 = - ^ + ^~?| «(A,0, ( 4 - 26 ) 
A-z(C- ^) 

where a (A , £) is given by Eq. (4.20). This completes the derivation of the boundary 
S matrix. 

S z = states 

We have already succeeded to determine the boundary S matrix. Nevertheless, 
a good check on this result and on the general formalism is provided by analyzing 
the S z = states, of which there are two. In particular, we consider the S z = 
state consisting of two holes in the Fermi sea, and also one 2-string. For £± — > oo, 
this is the spin-singlet (S = S z = 0) state shown in Fig. 6. * The position A of 



* The other S z = state is the one which for £± — > oo is one of the spin triplet 
(S = 1) states. For £± ^ oo, it is not clear how to identify this state in terms of 
the Bethe Ansatz solution, and we do not consider it further. 
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the center of the 2-string is not given by the simple expression (2.9). For example, 
for the special case £± = oo, the center position is 



(Ai) 2 + (A 2 )2 . (4.27) 



The general case £± 7^ 00 is discussed in Ref. 20. 

For the S z = states, the quantization condition (4.12) leads to a 2 x 2 matrix 
equation. The two eigenvalues of this matrix are pure phases. Since the matrix 
elements of R(X) and K(\,£±) are known, these eigenvalues can be computed 
explicitly. Let expi<E>(°) be the eigenvalue which for £± — > 00 corresponds to the 
spin-singlet (S = S z =0) state. The quantization condition implies 

e i2p(Ai)tf e i*<°> = 1 _ ( 4 .28) 

From Eq. (4.18) and the corresponding function r(A), we obtain the consistency 
condition 

e i(»(o)_»(D) = ei(Xi _ Aq) ei(Xi + Aq) ^ (4 29) 

where 

ei(A) = ^|, (4.30) 
A _ 2 

is defined in (4.17), the hole rapidities Ai and A 2 are arbitrary, and A is the 
corresponding rapidity of the center of the 2-string. 

For the case £± = 00, this relation is satisfied by virtue of the algebraic identity 

ei f J(Ai - A 2 )^ ei f l(Xi + A 2 )l = ei(Ai - A ) e^Ai + A ) . (4.31) 



which is true for arbitrary values of Ai and A 2 , where Ao is given by (4.27). We 
have explicitly verified the formula (4.29) also for the case £_ = 00, 7^ 00, and 
presumably it is true in general. This equality provides a nontrivial consistency 
check of the bulk and boundary S matrices and of the general formalism. 

4-3 Further remarks 

The boundary S matrix K(\,£) given by Eqs. (4.15), (4.20), (4.26) satisfies 
boundary unitarity and boundary cross- unitarity 18 , as well as the boundary Yang- 
Baxter equation 4 ' 18 ' 41 ' 42 

J R 12 (A-A / )^i(A,0 J R2i(A+A / )^ 2 (A / ,0 = ^2(A , ,0 J Ri 2 (A+A / )^i(A,0 J R2i(A-A / ). 

(4.32) 
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Since the bulk S matrix coincides with that of the sine-Gordon model with 
f3 2 — > 8/T, we expect that the boundary S matrix if(A,£) should coincide with 
the boundary £ matrix of Ghoshal and Zamolodchikov 18 for the boundary sine- 
Gordon model with (3 2 — > 8tv and with "fixed" boundary conditions. (For "fixed" 
boundary conditions, the field theory and hence the boundary S matrix are U(l) 
invariant.) We have verified that the two boundary S matrices indeed coincide, up 
to a rapidity-independent scalar factor, and with some redefinitions of variables. 
The bootstrap result of Ghoshal and Zamolodchikov for the boundary sine-Gordon 
model with "fixed" boundary conditions has been verified using the physical Bethe 
Ansatz approach by Fendley and Saleur 21 . Very recently, the boundary S matrix 
for the anisotropic spin 1/2 chain has been calculated by Jimbo, et al. 43 using the 
vertex operator approach. In the isotropic limit, their result coincides with ours. 

We have seen that the analysis of the S z = states for the open spin chain 
differs significantly from that of the closed spin chain. Indeed, for the open chain, the 
position of the center of the 2-string is a complicated function of the hole rapidities 
Ai and A2 (as well as the boundary parameters £±); while for the closed chain, the 
center of the string is located midway between the two holes. Naively, one might 
worry that this leads to a breakdown of factorization. However, we have seen that 
factorization is maintained by virtue of certain nontrivial identities. We expect that 
a similar situation holds for the closed spin chain with four or more holes. 
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Figure Captions 

Fig. 1: Ground state of spin 1/2 chain, with iV = 30. Diamonds denote real 
roots of the Bethe Ansatz equations. 

Fig. 2: S = S z = 1 excited state. Open circles denote holes. (The scale here 
differs from the one in Fig. 1.) 

Fig. 3: S = S z = excited state. The 2-string (denoted by X's) has its center 
at (Ai + A 2 )/2. 

Fig. 4: Two-dimensional vertex model. Solid and dashed lines correspond to 
spin 1/2 and spin 1, respectively. The vertex formed by lines corresponding to 
spins si and S2 has weight R( Sl ' S2 \\). (Periodic boundary conditions should 
be imposed in both horizontal and vertical directions.) 

Fig. 5: Ground state of alternating spin 1/2 - spin 1 chain. Diamonds denote 
real roots (1-strings) and X's denote complex roots (2-strings) of Bethe Ansatz 
equations. 

Fig. 6: S = S z = excited state. The center of the 2-string is given by 
Eq. (4.27). 
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